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We report the observation of vortex nucleation in a rotating optical lattice. A ^^Rb Bose-Einstein 
condensate was loaded into a static two-dimensional lattice and the rotation frequency of the lattice 
was then increased from zero. We studied how vortex nucleation depended on optical lattice depth 
and rotation frequency. For deep lattices above the chemical potential of the condensate we observed 
a linear dependence of the number of vortices created with the rotation frequency, even below the 
thermodynamic critical frequency required for vortex nucleation. At these lattice depths the system 
formed an array of Josephson-coupled condensates. The effective magnetic field produced by rotation 
introduced characteristic relative phases between neighbouring condensates, such that vortices were 
observed upon ramping down the lattice depth and recombining the condensates. 

PACS numbers: 03.75.Lm, 67.85.Hj, 74.50.+r, 74.81.Fa 



Ultracold quantum gases in rotating optical lattices 
find themselves at the intersection of two fields which 
have generated an impressive body of experimental and 
theoretical work. The versatile and clean potentials of- 
fered by optical lattices have proved to be an incredi- 
bly adept system for exploring a wide range of funda- 
mental problems in condensed matter physics such as 
the Mott-insulator transition |T], Anderson localization 
[2113] and Tonks- Girardeau gases [4]. Similarly rotating 
quantum gases have provided an array of striking results, 
from the nucleation of vortex lattices in Bose gases [3 
to rotating Fermi gases at the BCS-BEC crossover [7|. 
The close analogy between the physics of rapidly rotating 
neutral atoms and electrons under a magnetic field has 
led to considerable interest in the possibility of achiev- 
ing strongly-correlated quantum Hall states in a rapidly 
rotating atomic gas [H [9] . 

Rotating optical lattices for cold atoms have gener- 
ated a large amount of interest in their own right, which 
can broadly be divided into three main areas: (i) Weak 
rotating lattices have been used to pin vortices [10 and 
theoretical work has shown that rich vortex lattices struc- 
tures, including lattices of doubly quantized vortices, are 
predicted to emerge [Til HI]; (ii) Stronger optical lat- 
tice potentials with a large number of atoms per site 
(100-1000) realize the physics of Josephson junction ar- 
rays ( J JA) under magnetic fields [131 HI] ; (iii) Arguably 
the most interesting regime is for a dilute gas in the 
tight binding regime in a rotating lattice where fractional 
quantum Hall physics is predicted to occur p!5HT7] . In 
all three regimes rich structures emerge that depend on 
the relative density of vortices and lattice sites. 

We report on the first experiments with a rotating op- 
tical lattice in the 'deep' lattice regime where lattices 
depths were reached such that a 2D array of weakly linked 
condensates was created, forming a Josephson junction 
array. Under rotation this system realizes the uniformly 



frustrated J JA [131 HH] • We have investigated the nucle- 
ation of vortices in a rotating lattice, starting from a non- 
rotating condensate loaded into a static lattice, and then 
increasing the rotation rate of the lattice. There has been 
one previous experiment employing a rotating optical lat- 
tice by Tung et al. [lOj at JILA, where it was reported 
that heating due to mechanical instabilities and aber- 
rations limited the rotating lattice to depths less than 
30% of the condensate's chemical potential. We were 
not limited to this weak lattice regime and our system is 
markedly different to the JILA experiment where a vor- 
tex lattice was created by an elegant evaporative spin-up 
technique [6] before a weak co-rotating optical lattice was 
imposed. 

The rotating lattice used in the experiments described 
in this paper was generated by dual-axis acousto-optic 
deflectors and a novel optical system, described in detail 
in [19 . The two-dimensional optical lattice was formed in 
the focal plane of a custom-made multi-element lens with 
a numerical aperture of 0.27. Generating the rotation 
of the lattice acousto-optically gave intrinsically smooth 
rotation and allowed the rotation rate to be varied during 
an experimental run. 

The Coriolis force on neutral atoms in the rotating 
frame has the same form as the Lorentz force experi- 
enced by charged particles in a magnetic field. Rotation 
at angular frequency ft gives rise to the effective mag- 
netic vector potential A = 11 x r, and the Peierls phase 
gained by an atom tunnelling from site i to site j is 
Aij = (m/h) J^^ A(r') • dr', where m is the mass of the 
atom. The sum of the link phases around a lattice pla- 
quette of side d obeys the property Zlpiaquette = ^^Z' 
where / = 2mQ(P /h is the mean number of vortices per 
lattice plaquette. 

The experiments started with a ^^Rb condensate con- 
taining 2 X 10^ atoms with no visible thermal compo- 
nent, held in a trap formed by an axial symmetric loffe- 
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Pritchard trap and a single red-detuned light sheet in 
the radial plane increasing the axial trapping frequency 
to give {ujr^^z} — 27r{20.1, 53.0} Hz. The residual ra- 
dial anisotropy of the combined magnetic and optical 
trap was uoy/ujx = 1.008 ± 0.003. The dipole sheet trap 
was creating using broadband light with a spectral width 
AA ~ 3 nm centred on 865 nm in order that the light 
had a short coherence length. This ensured there were 
no interference effects, e.g. from multiple reflections at 
the vacuum cell wall, adding noise to the trapping po- 
tential. With narrow-band light it was observed such 
perturbations acted to spin down the cloud and caused 
irregularities in the filling of the lattice. 

The 2D optical lattice was formed in the radial plane 
of the trapped condensate by four circularly polarized 
beams at A = 830 nm intersecting in the focal plane of a 
high N.A. objective lens, achieving a lattice constant of 
2 /im. One pair of beams was detuned by 10 MHz to the 
other pair ensuring interference between the orthogonal 
optical lattices did not affect the atoms. The optical 
lattice was initially static as it was ramped to its final 
depth Vb in the range 100 Hz < < 4000 Hz [20 . The 
radial Thomas- Fermi radius of the condensate was 16 /im, 
resulting in the filling of ~ 16 lattice sites across the 
diameter of the condensate and around 200 sites in all. 
The lattice enhanced interaction energy at the central 
sites rose above the chemical potential of the unperturbed 
condensate (/i ~ 500 Hz) meaning an array of individual 
condensates was not formed until Vb ~ 1000 Hz. The 
central wells contained around 1500 atoms. 

After the condensate had been loaded into the static 
optical lattice was increased linearly from zero to its 
final value in 320 ms, followed by a further 100 ms of 
rotation at the final value. While the lattice was still ro- 
tating the lattice depth was then ramped down to zero 
in 12 ms, allowing the condensates to merge and con- 
verting phase differences around plaquettes to vortices. 
Previously Scherer et at [21] have investigated vortex for- 
mation by merging three independent condensates, and 
Schweikhard et al ^2] have observed the thermal activa- 
tion of vortices on a two-dimensional lattice of Josephson- 
coupled BECs, however our experiment takes place in 
the rotating frame. Rotation gives rise to an additional 
phase 27r/ around a plaquette of the optical lattice, sim- 
ulating the effect of a magnetic field on charged particles. 
After the optical lattice was ramped down the cloud was 
then immediately released from the magnetic and optical 
traps and destructively imaged after 20 ms time- of- flight 
expansion. We also repeated the sequence but holding 
the cloud after the lattice was ramped down for an ad- 
ditional 500 ms, to allow any vortex- ant ivort ex pairs to 
annihilate (such pairs could be created in a static lat- 
tice [22 ). This generally improved visibility of vortices 
in the harmonic trap and we observed the same number 
of vortices in the system within experimental error. The 
uncertainty in the number of vortices observed for a given 
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FIG. 1: (a) Number of vortices nucleated as a function of 
optical lattice rotation frequency Q for two lattice depths, 
Vq = 190 Hz and Vq = 1200 Hz. The error bars denote 
the standard deviation in vortex number for three experi- 
mental images. Inset: an example of a vortex lattice created 
by the rotating optical lattice of depth Vb = 1200 Hz at Q./ 
27T = 10 Hz. (b) The minimum rotation frequency needed 
to nucleate a single vortex, Qmin, as a function of the lattice 
depth Vb (1 Hz was the minimum frequency used). 



ft arose from difficulties in determining the presence of 
vortices near the edge of a condensate. 

Figure[TJa) shows the number of vortices, Ny, observed 
as a function of the optical lattice rotation frequency, 1], 
for lattice depths of Vb = 190 Hz and Vo = 1200 Hz. The 
behaviour of Ny{Q) has a characteristic shape which de- 
pends on the exact vortex nucleation process. The first 
BEG stirring experiments used stirring beams of a similar 
size to the condensate and observed that no vortices were 
nucleated below ft ^ O.To;^, where there is excitation of 
the / = 2 surface mode [23 . It was later shown that suf- 
ficiently small stirring beams nucleate vortices without 
discrete resonances of surface modes [24 , and the mini- 
mum frequency for vortex nucleation, ^Imin, matched the 
frequency at which a single vortex at the center of a con- 
densate is energetically stable, ftc [25^- For om experi- 
mental conditions 17c/27r ~ 4 Hz. For Vb = 190 Hz, far 
below the ~ 1000 Hz needed to split the condensate into 
weakly linked islands. Fig. [TJa) shows a similar trend 
to that observed in [24], i.e. Ny did not exhibit sharp 
resonances and l^min ^ ^c- We checked this behaviour 
was an equilibrium property by repeating the experiment 
with rotation times up to 1 s, and observed the same re- 
sponse. We infer that a weak (Vb ^ /i) rotating lattice 
stirs up a condensate in a non-resonant fashion. For a 
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FIG. 2: Number of vortices nucleated as a function of Q/ 
ujr. The solid line is the number of vortices expected for a 
cloud of fixed radius, i^o = 16 /xm, while the dashed line in 
the number of vortices expected for a centrifugally distorted 
cloud of radius R = Ro[l - (Q/cj^)^]-^/^^. We observed the 
optical lattice prevented the cloud from reaching R for a given 
allowing rotation above the critical frequency Q > uJr to 
be achieved. 

rotating lattice of depth Vq = 1200 Hz, however, Ny 
displayed a distinctly different trend, with the number 
of vortices observed following a linear dependence on Q 
even below Qc This change in behaviour is indicative 
of entering the regime where the system forms an array 
of Josephson-coupled condensates, rather than simply a 
condensate perturbed by a weak rotating lattice as for 
Vb = 190 Hz. Figure [TJb) highlights the transition be- 
tween the two regimes, showing the dependence of l^min 
on Vq. l^min f^Hs bclow Qc as Vb rises above /i = 500 Hz, 
and falls to 1 Hz (the minimum rotation frequency used) 
for Vb > 1000 Hz, corresponding to the lattice depth 
at which condensates are expected to be well localized 
on sites and can only communicate by tunneling. The 
rotation induced phase differences between neighbouring 
condensates are converted to vortices when the lattice 
is ramped down in the rotating frame. This nucleation 
mechanism creates vortices near the centre of the con- 
densate even when Q < Qc- Vortices do not have to spin 
in from the edge of the cloud, which is what we observed 
for stirring with a weak lattice potential. 

The measurements of the number of vortices Ny was 
extended to higher rotation frequencies Q as shown in 
Figure [2] The ratio Q/uJr is used where uJr is the ra- 
dial trapping frequency taking into account the slight en- 
hancement by the Gaussian envelope of the lattice beams: 
for Vb = 600 Hz, cj^/27r = 21.5 Hz and for Vb = 1100 Hz, 
co'r/27r = 22.6 Hz. The density of vortices for a given 
whether a lattice is present or not, is Uy = mQ/nh. The 
total number of vortices is then Ny = mQR^ /h. For a ro- 
tating BEG in equilibrium at Q the centrifugal distortion 
of the Thomas-Fermi radius is R = Ro[l - {ft/ujr)'^]-^^^^ 
[26] . The predicted number of vortices for a cloud dis- 
playing this behaviour is denoted by the dashed line in 




FIG. 3: Images after release from the rotating optical lattice 
showing the pattern of vortices: (a) Gondensate before rota- 
tion; (b) Q/27V = 4 Hz; (c) Q/27V = 7 Hz; (d) Q/27V = 17 Hz. 

Fig. |2] The solid line presumes the radius stays constant 
at i?o = 16 /im. For Vb = 1100 Hz, A/"^ rises linearly with 

whereas for the weaker Vo = 600 Hz lattice Ny shows 
a faster, nonlinear dependence on Q. We attribute this 
behaviour to the deeper lattice suppressing the effect of 
centrifugal distortion. For the deeper lattice atoms may 
only redistribute themselves by tunneling between sites. 
While the tunneling parameter (the Josephson-coupling 
energy [13]) is around 2 kHz at central lattice sites this 
drops to < 50 Hz at peripheral lattice sites due to the 
drop in atom number, giving tunneling times at the edge 
of the lattice comparable to the experimental duration. 
Bloch oscillations may also play a role in inhibiting the 
atoms from spreading out in the lattice [27 . 

As a result of the optical lattice inhibiting the spread- 
ing out of the cloud we were able to rotate the cloud 
above the critical frequency fl > oUr- As seen in Fig. [2] 
the number of vortices dropped for rotation above cj^, due 
to increased heating of the condensate. However vortices 
could be observed up to 1] = 1.15 cj^ for Vo = 600 Hz 
before heating and reduction of vortex visibility meant it 
was no longer possible to observe vortices. The heating 
dit Q ^ UJr was worse for the deeper Vb = 1100 Hz lattice 
and vortices could not be observed above Q = 1.02cj^. A 
cloud in a harmonic trap without the presence of an op- 
tical lattice should be expelled for Q > Ur though Bretin 
et al. were able to achieve Q = 1.06 uJr using additional 
quartic confinement [28 . 

The structure of the vortex patterns formed in a 
bosonic JJA under rotation are heavily influenced by the 
presence of the lattice. For values of the frustration pa- 
rameter (the mean number of vortices per lattice plaque- 
tte) / = 2mVtd? /h = p/q (where p and q are integers) the 
ground state vortex patterns have been calculated [13], 
and have a unit cell of size q x q sites. Between these 
rational fraction values of / intermediate patterns are 
formed. Figure |3] shows examples of the vortex patterns 
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FIG. 4: Condensate fraction Nq/N as a function of lattice 
depth Vb after ramping up the rotation rate to 15 Hz in 40 ms 
followed by 35 ms of rotation at 15 Hz. 



we observe. We do not reproducibly observe the pre- 
dicted q X q unit cell structure as we scan Q and hence 
/. There are several possible reasons for this: (i) for 
d = 2 jam and uJr/^TT = 20.1 Hz the maximum value of 
/ we can reach is / ~ 1/5. The finite size of the system 
then becomes an issue, as it will be difficult to realize 
the periodic structure for q = 5 when a maximum of 200 
sites are filled by the condensate; (ii) the energy land- 
scape of cold atoms in a rotating lattice is known to be 
very glassy, that is closely spaced minima in energy can 
be separated by large barriers, making it hard to reach 
the true ground state Lack of equilibration on ex- 
perimental timescales might then account for the absence 
of qx q vortex structures; (iii) any intrinsic heating from 
the rotating optical lattice will make it difficult to reach 
the ground state. 

The issue of finite size can be addressed by increasing 
Uz or using a larger condensate, enabling more lattice 
sites to be filled. The values of / achievable can be in- 
creased by implementing tighter radial trapping to allow 
higher Q to be reached or increasing the lattice constant 
d. The / = 1/2 case is particularly interesting due to 
competition between an Ising-type and BKT type phase 
transition [14]. Vortex patterns in a bosonic JJA under 
rotation for < / < 1/2 will be explored in future work. 

An important question regarding the viability of fu- 
ture rotating optical lattice experiments is the issue of 
heating. Figure [4] displays the fraction of atoms in the 
condensate, Nq/N, after rotation at Q/27r = 15 Hz {Q/ 
27T was increased from to 15 Hz in 40 ms followed by 
35 ms of rotation at 15 Hz) for different lattice depths, 
showing an approximately linear dependence. We found 
at Q/27r = 15 Hz, Vq = 1100 Hz we could rotate for 
around 300 ms before heating significantly reduced the 
condensate fraction, with longer rotation times achiev- 
able for smaller or or both. The main sources of 
heating were fiuctuations in the lattice depth during ro- 
tation (at the ±1 % level) and residual imperfections in 
the lattice upon rotation. All the data in this paper were 
taken without using any additional evaporation during or 
after the rotation process (unlike previous stirring experi- 
ments [5l|23]). The implementation of an additional cool- 



ing scheme in the optical lattice such as that described 
by Griessner et al. [30 could enable sufficiently low tem- 
peratures to be achieved to realise quantum Hall physics 
in future experiments with the rotating lattice p!5HT7] . 

In conclusion, the use of acousto-optic defiection to 
generate smooth rotation of an optical lattice at a well- 
defined frequency has allowed us to explore a new method 
of vortex nucleation. For lattice depths Vq < fi we found 
the weak lattice acted as a stirring mechanism. For 
deeper lattices the system was split into well-localized 
condensates at each site and a linear dependence of vor- 
tex number on Q was observed for rotation frequencies 
below the thermodynamic critical frequency for vortex 
nucleation. This implied vortices were created locally at 
a lattice plaquette. 

This research is supported by EPSRC, QIP IRC 
(GR/S82176/01) and ESF. 
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